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Abstract 

The regular solutions for the Ginzburg-Landau (-Nielsen-Olesen) Abelian 
gauge model are studied numerically. We consider the static isolated cylindri- 
cally symmetric configurations. The well known (Abrikosov) vortices, which 
present a particular example of such solutions, play an important role in the 
theory of type II superconductors and in the models of structure formation 
in the early universe. We find new regular static isolated cylindrically sym- 
metric solutions which we call the type B and the flux tube solutions. In 
contrast to the pure vortex configurations which have finite energy, the new 
regular solutions possess a finite Gibbs free energy. The flux tubes appear to 
be energetically the most preferable configurations in the interval of external 
magnetic fields between the thermodynamic critical value H c and the up- 
per critical field H C2 , while the pure vortex dominate only between the lower 
critical field H Cl and H c . Our conclusion is thus that type B and flux tube 
solutions are important new elements necessary for the correct understanding 
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of a transition from the vortex state to the completely normal state. 
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1. INTRODUCTION 



The Ginzburg-Landau theory of superconductivity |T| is mathematically equivalent to 
the Abelian theory of coupled gauge and Higgs fields 0. The existence of vortex (string- 
like) solutions in it was predicted by Abrikosov in 1952 (and published 5 years later ||) in 
the context of the phenomenological model of superconductors and was discussed in in 
the framework of the dual string approach (see also a recent generalization to the case of 
nontrivial helicity in M). The aim of this work is to clarify, with the help of the careful 
numerical analysis, the properties of such solutions and to construct new solutions. We 
confine ourselves to the case of cylindrical symmetry, thus considering an isolated vortex and 
related field configurations. In the literature devoted to this subject (see, e.g., [pt^|,^HT0| ) 
main attention was paid to approximate solutions and qualitative methods, but there were 
very few attempts to study exact solutions by numerical methods. Moreover, as recently 



was shown in ||11|| , the old qualitative results may contain mistakes and are incomplete. 
Approximate results are normally confined to the domains of very small (A 1) or extremely 
large (log A ^> 1) values of the coupling constant A (or the characteristic Ginzburg-Landau 
parameter k = see Appendix). As for numerical studies, one can mention the earlier 
reports 112-1^1 and more recent (variational) analysis in [H]. In our paper we present 
new numerical results for the isolated regular structures in the Ginzburg-Landau (-Nielsen- 
Olesen) model. Although here we confine ourselves to the case of Abelian gauge field, we 



will consider elsewhere the non- Abelian generalizations (cf. previous approaches in [iq,[L9|0. 
Our results demonstrate a rich structure of the space of exact solutions for the classical 
Ginzburg-Landau model. Besides the well known vortex solutions, for which we compute a 
variety of parameters in an interval of values of A close to the critical value A = | (k — \J\), 
we describe some new exact regular solutions. In our opinion, of particular interest are the 
configurations which we call the flux-tube and the oscillating solutions below. The former 
appear to comprise a new nontrivial structure in ideal type II superconductors which is 
important for understanding the transition between a pure vortex state and a normal state. 
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As it is well known, the vortex (or the mixed) state of a type II superconductor exists 
between the critical magnetic fields H Cl and H C2 . Their values, as is usually claimed in 
the literature []3|J5| |8|,[T0f] , can be determined from the analysis of isolated vortex solutions. 



However, from our results we conclude that such a claim is correct only partly. Indeed, 
H C1 is determined, for each value of A = k 2 , by a relevant isolated vortex solution from 
a comparison of the Gibbs free energy of the vortex configuration with that of the pure 
Meissner configuration. But the upper critical magnetic field H C2 appears to be completely 
unrelated to the vortex configurations (cf. ||). Instead, one can only determine H C2 from 
a similar comparison of the Gibbs free energy for a normal state with that of a different 
type of exact regular solutions of the Ginzburg-Landau equations. We call these the type 
B solutions. Actually, approximate type B solutions are usually described in the literature 
PJ5|-[T0| in the framework of analysis of the linearized Ginzburg-Landau equations. We 



construct numerically exact self-consistent type B solutions for all possible values of external 
magnetic field. Although interesting enough themselves as mathematical structures, the type 
B solutions have a limited physical value as compared to the new flux tube solutions. The 
latter configurations correctly describe the transition from a pure vortex to a normal state, 
as will be demonstrated below. 

It seems worthwhile to notice that the Ginzburg-Landau equations describe common 
extremals for two different (action type) functionals: for the usual energy and for the Gibbs 
free energy integrals. The principal difference of the vortex and the new solutions is that the 
former are the finite energy configurations, while the latter are the finite Gibbs free energy 
configurations. Previously, attention in the literature was paid only to the finite energy 
regular solutions of the Ginzburg-Landau equations. In the context of cylindrical symmetry, 
these are the famous Abrikosov(-Nielsen-Olesen) vortices. It is our aim to demonstrate the 
existence and the physical relevance of two large families of finite Gibbs free energy regular 
solutions which are described in our paper as the flux tube (or type A) and the type B 
solutions. In our opinion, the correct understanding of the mixed state of type II supercon- 
ductor can only be achieved after taking into account these new solutions. In particular, 



it turns out that the class of the flux tube solutions is divided into an infinite number of 
families labeled by a number of nodes n for the scalar field configuration. Each family exists 
on a finite interval of magnetic field (definition of limiting points see below). In this way, 
one finds a certain fine structure of the mixed state for an ideal type II superconductor. 

The plan of the paper is as follows: Section 2 contains a general introduction to the 
model, we discuss the two energy functionals and formulate general regularity conditions 
at the origin. A brief account of numerical analysis of the pure vortex solutions is given in 
Section 3. Sections 4 and 5 contain the description of new solutions, the type B and the 
flux tubes, respectively. In Section 6 we perform the linearization analysis of the Ginzburg- 
Landau equations, while Section 7 is devoted to the regular oscillating solutions. The latter 
appear to be certain unstable "relatives" of the vortices. Finally, Section 8 contains a 
short discussion and a summary of the results obtained. In the Appendix we compare our 
notations and conventions with the old ones used in the literature. A general remark is 
necessary for the tables and figures: We find it convenient to put all the available numerical 
data into the separate addendum EO] , available in electronic form. In the present paper, the 



Tables contain only some selected reasonable minimum of data, while the Figures present 
additional information. 

2. ENERGY FUNCTIONALS AND REGULARITY CONDITIONS 

A. Nielsen-Olesen Lagrangian 

The Lagrangian of the Abelian ("Nielsen-Olesen" 0]) gauge model describing interacting 
electromagnetic F = dA and complex scalar field $ reads 



with the potential 



£no = ~^{F A *F + D$ A *£>$) - V(|$|)77, (2.1) 



V(m) = ^(m 2 -^) , (2-2) 



where the overbar denotes complex conjugation and D = d + iA, and star * is the Hodge 
dual operator, n :=*1 being the volume form. 

Noticing that p has a dimension of inverse length, we can introduce for the cylindrical 
system (p, 9, z) a new (dimensionless) radial coordinate r via p = —r. We are looking for 
static configurations, and use the following cylindrically symmetric ansatz (cf. 0) 

A = f(r)d9, *=JL<p(r), (2.3) 

where /, ip are two real functions. 

The magnetic field 1-form, defined as the (three-dimensional) Hodge dual *F, has only 
one component in z-direction, *F = Hdz, where the latter is given by the expression 

A r dr 

We will denote the dimensionless magnetic field by h := HX/p 2 = f'/r. 

Magnetic field is conveniently characterized by the flux it produces through a two- 
dimensional surface. In a cylindrically symmetric case the total flux over a surface orthogonal 
to the z-axis is 

F — [ pdpd9 H = 2n r dr%- = 2tt[/(oo) - /(0)1. (2.5) 
J Jo dr 



B. Energy (line density) functional and equations of motion 



For the Lagrangian (2.1) one can immediately write the energy per unit length in the 
form 



£ 



oo ] 

dr r— 
X Jo 2 



\r dr J \dr J r 2 ^ 2^ ' 



(2.6) 



Notice that it is not necessary to include a phase factor e md for the scalar field in the ansatz 
([2.3p, as it is done in a number of different approaches. The field <p is always defined up to 
a gauge transformation and we find it more convenient to work in the gauge (|2.3| ). 
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Now, it is straightforward to see that the equations of motion of the Nielsen-Olesen model 
read: 

r 2 f" - rf = r 2 p 2 f, (2.7) 
rV" + np' = <p(f 2 + Ar V - 1)). (2.8) 

The value of the constant A plays an important role. If it equals to the critical value A = ~, 
the Nielsen-Olesen equations are consequences of the first order (Bogomolny) system 

V + ^ 2 "1)=0, (2.9) 

<ff+-f<p = 0. (2.10) 
r 

In the theory of superconductivity A = \ separates two phases: for A > | (resp., A < |) one 
has a type II (resp., type I) superconductor. The critical subcase was extensively studied in 
the literature [pT H23|JTT|1 . The general noncritical case for A ^ | is much less investigated. 



Using ( |2.4| ), we can write the equation fl2.7| ) in the form 



h' = -fp 2 , (2.11) 
r 



and hence the system (|2.7|)-([2.8|) can be transformed into 



h" + ~ti = V 2 + 2h'^-, (2.12) 



1 , 



+ = v \- i (h'Y + x{^-i)y (2.i3) 

explicitly for the coupled magnetic field h and scalar field (p variables. This system is however 
not particularly useful for numerical investigation because of explicit l/<p terms. 



C. Gibbs free energy functional 

The line energy density functional ( |2.6| ) is minimal (£ m in = 0) for the Meissner state, i.e., 
when ip = 1 (superconducting order in all points of a sample) and h = (no magnetic field 
inside a superconductor). This corresponds to a trivial solution of (|2.7| )- (|2.8| ), / = 0, (p — 1. 

7 



The normal state is described by another simple solution, f = C + \h r 2 , ip = 0, (C and h 
are constants), which always (with or without magnetic field ho) has formally infinite energy 
S. However physically, of interest is a difference of energies, not an energy itself. 

In particular, let us turn our attention to the Gibbs free energy (per unit length) which 
is defined by 



g = s 



I 



d 2 x (HH 



ext / 



(2.14) 



We assume that both magnetic fields, the internal one H and an external H ext are directed 
along the z-axis, and use the dimensionless values of fields defined, in accordance with ( |2.4| ), 
by H = (/j 1 2 /X)h and H ext = (fi 2 / X)h cxt . The Gibbs free energy of the normal state with a 
magnetic field h = h ext is given by the integral 



Q nh = £[ip = 0,h = h ] 



2tt// 2 



drr hoh e 



2-nn 2 f°° , 1 
— - — / drr - 
X Jo 2 



(2.15) 



For a sample without boundary, this is an infinite constant, while for a cylinder of a radius 
R this is a finite positive (negative) constant irR 2 fi 2 (l/2 — h\fX)j2 for a magnetic field 
ho < h c (h > h c ), and zero for h = h c :- 



'A/2. This observation underlies the physical 
interpretation of h c as a critical (so called thermodynamic) value of the magnetic field which 
distinguishes normal and superconducting phases. 
Now, let us consider the difference, 



A£ := Q — Q n h 
2tth 2 i-oo 



drr — 
X Jo 2 



'ld£ 
r dr 



2 / , \ 2 



1 



;2.i6) 



It is very important to notice that the Gibbs functional ( |2.16| ) has the same equations for 
extremals as the energy functional ( |2.6| ), namely ( p.7| )-( p^8| ). However, unlike the strictly 
positive (^.6|), the functional ( |2.16| ) can have any sign. 



D. Regularity at the symmetry axis 

Looking for solutions which are regular at the origin, we substitute the series expansions 



S 



/ = £a fc r fe , v = Y,hr k , (2.17) 

k=0 k=0 

into (|2.7Q - (|2T8|) . We then find two types of conditions: 

(A) Potential / is non-zero while scalar field <p vanishes at the origin, 

/ = N + or 2 + m „ + 1 f r™» + 0(r 2N n (2.18) 

^ 6r "( 1 + WTT)( a -^) r2 + 0(r4) )' (2 - 19) 

where N = ±1, ±2, ... is a nonzero integer and parameters a, b are arbitrary. 

(B) Potential / vanishes while scalar field p is nontrivial at the origin, 

/ = ar 2 (l + h) 2 r 2 + 0(r 4 )) , (2.20) 
V = b(l + ^(6 2 - l)r 2 + 0(r 4 )^) , (2.21) 

with some parameters a, b. 

In both cases, parameter a determines the value of the magnetic field at the origin, 
h(0) = 2a. When 6 = 0, both cases reduce to the solution which describes a normal 
superconductor if = filled by the constant homogeneous magnetic field h(r) = h(0) = 2a. 

3. VORTEX SOLUTIONS 

Vortex solutions of the Ginzburg-Landau equations (|2.7|)-(|2TE|) are distinguished among 
others by the special conditions at infinity which read: 

/(OUoo— 0, (3.1) 
^(tOUoo— >1. (3.2) 



When combined with the type A regularity conditions at the origin (|2.18|) -( p.l9[) , these 



asymptotic conditions define uniquely parameters a, b for any value of A and N. The meaning 
of the constant is clear: this is the value of the magnetic flux described by such solutions. 
Indeed, we substitute /(0) = N and /(oo) = into ( |2.5| ) to obtain F = —2ttN. This is the 
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well known flux quantization result for superconductors. Notice that the asymptotics ( |3.1| )- 
( p.2| ) cannot be realized for any solution with type B regularity conditions ( |2.20| )- ( pT2T| ) , see 
the discussion in the next section. 

The results of the numerical integration are given in Figure [I] and Table [I], where left 
and right columns describe solutions for N = 1 and N = 2, respectively. For convenience of 
comparison with the old calculations [|l2Hl6|l we write the coupling constant as A = k, 2 [see 
Appendix for definitions, in particular QA.5 )]. 



Energy (|2.6|) is always positive on all vortex configurations. It decreases with growing 
K and increases when k — > 0. A useful physical variable, which allows to find necessary 
conditions for the existence of vortices, is the Gibbs free energy (line density) ( |2.14|) . This 
quantity is zero for the Meissner state (h = 0, </? = 1), while for a vortex state it is always 
positive when an external magnetic field H ext is oriented oppositely to the magnetic field H 
inside a sample (hence such vortices are ruled out), and it can become negative for a certain 
value of external field with the same orientation as in the sample. In the latter case, for a 
constant external magnetic field H ext = (0, 0, H ext ) one finds 



Q = £ - 2nH { 



ext. 



drrh 



o 



£ - 2tt \N\ H ext . (3.3) 



Hence the lower critical value of an external magnetic field H cxt = H Cl , determined by the 
condition Q = 0, is equal to 

where in the last equality we switched to the old conventions summarized in Appendix. 

It seems worthwhile to make the following short remark: The analysis of exact vortex 
configurations shows that the actual (computed for vortices) magnetic field penetration 
length for different values of k is always greater than the domain of an essential change of the 
scalar field, even for small k. In fact, this observation was also confirmed in earlier numerical 
studies [|T3|,|T^]. This casts doubts on a possible qualitative and quantitative understanding 
(cf., [H-01]9) °f these solutions on the basis of so called penetration lengths 8 and £ (see the 
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Appendix) for the magnetic and the scalar fields, and stresses the need of exact numerical 
investigations (see details in addendum 0] where the numerous vortex solutions for a wide 
range of k are described). 



4. REGULAR SOLUTIONS FOR TYPE B CONDITIONS 



The type B conditions ( 2.20f) - (|2.21|) , like the type A conditions (|2.18|) - (|2.19 ), guarantee 



a regular behavior of solutions at the origin. (It seems worthwhile to notice that type 
B conditions are not a particular case of type A for N = 0). It is straightforward to 
analyze a qualitative behavior of such solutions for finite values of r. For example, equation 
( [2.7|) immediately yields that f(r) is a monotonous increasing for a > (decreasing for 
a < 0) function. Indeed, if we assume an extremum for a finite r = r , one finds from ( |2.7| ) 
f"{ r o) — y? 2 /( r o) which means that such an extremum is a minimum for positive /(r ) and a 
maximum for negative /(r ). For a function starting from a zero ( |2.20| ) both possibilities are 



excluded and thus f(r) is a monotonous function. In particular, this means that f(r) — > oo 
for r — > oo, and one thus concludes that no finite energy solutions for type B regularity 
conditions exist: the functional ( |2.6|) is infinite. 

Let us however look at Q2.7D - ([2.8D as the equations for extremals of the Gibbs free energy 
functional (|2.16| ). It is immediately clear that there exist finite Gibbs free energy regular 



type B solutions, provided they satisfy at infinity 



h(r) = -f'(r) 
r 



— ho, (4.1) 

• oo 

<^(r)Uoo^0. (4.2) 



Notice that unlike in a vortex configuration (|3.2j ), the scalar field cannot approach at infinity 
any finite value except 0. 

The results of numerical integration are presented in Tables 0-fT| and Figure 0. A 
principal difference of these solutions (which we will call "type B" solutions for brevity) 
from the vortex solutions lies in the fact that the magnetic field is asymptotically constant, 
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and hence the flux integral defined by (|2.5| ) is formally infinite. However a reasonable 
replacement is provided by the quantity 



M 



AF = J pdpdd (H - H ) = 2n J drr(h-h ). (4.3) 



Defined formally as a difference of fluxes, this variable is usually interpreted as a magne- 
tization per unit volume [f7Ul4"|. Unlike the quantized flux for the vortices, M can have an 
arbitrary value. 

From Tables [Hj-pl we see, for the type B solutions without node, that when a — > A/2 
the external magnetic field is approaching h ext = h C2 = A, while M/{2tx) = AQ = 0. For 
A > 1/2, we find a higher external magnetic field for larger values of magnetic field at the 
center h(0) = 2a. When A = 1/2, for all solutions h ext = 1/2 (see JZIJ). Finally, for A < 1/2 
|2(| the external magnetic field is decreasing when h(0) = 2a grows (the reverse order as 
compared to the A > 1/2 case). 

Besides the simple type B solutions without nodes, there exist more nontrivial solutions 
with nodes. Both are displayed on Figure ^ while the configurations with one node are 
represented by the right columns in the Tables 0-|T|. 

For A > 1/2 all type B solutions without node have a negative Gibbs free energy AQ, 
for A = 1/2 always AQ = 0, and for A < 1/2 the Gibbs free energy is positive. One can find 
also negative AQ for type B solutions with one node when h c < A/3 (i.e., k > 3/y/2), see 
the case of k = 2.25 (Table pT[) . This will be clarified later in Section 6. 

The magnetization curves for the type B solutions are given in Figure [|. For any A = k 2 , 
M diverges at h c . 

In the present paper, we display the results for k = 1.0,2.25, the detailed tables and 
figures for k ranging from 0.5 to 5.0 are collected in the addendum pi . 



5. FLUX TUBE SOLUTIONS 



Let us look for the other solutions which yield finite values for the Gibbs functional 
( |2. 161 ). At the origin r = we take the type A regularity conditions ( |2.18| )-( p.l9| ), while at 
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infinity we consider the asymptotics (|4.1|) -( p~2j ). Such a combination of two type conditions 
at zero and at infinity suggests a possible physical interpretation of such solutions which 
appear as a result of a certain "gluing" of a vortex configuration at the origin with a type 
B solutions at large radial values. 

The results of numerical integration are given in Figures [||]| and Tables fVJ-[V|. As in the 
previous section, we present results only for the values of k = 1.0, 2.25, see |2(| for more 
solutions. 

Figure ^| explains why we call these solutions the flux tubes: There is a core where matter 
is in a state close to the normal one filled by the magnetic field (this is in fact a vortex), 
surrounded by a superconducting tube (almost completely free of a magnetic field). Outside 
such a tube the sample quickly reduces to a normal state with the external field penetrated 
in it. For the solutions with one node we have a "sandwich-like" structure: a tube of normal 
state between two superconducting tubes. 

Each family of flux tube solutions has two branches. E.g., for k — 1, one of these 
branches is characterized by the positive Gibbs free energy, and another has negative Gibbs 
free energy. However, for sufficiently large k both branches describe the negative Gibbs 
free energy configurations. We find it convenient to depict these branches in the form of 
magnetization curves, Figure £|. For A < |, two branches with positive Gibbs free energy 
exist. As in the case of type B solutions, M diverges at h c . 

It is interesting to notice that among the flux tube configurations there are solutions 
with quite unusual behavior of the magnetic field, which in the center is oppositely oriented 
with respect to the direction of external magnetic field. See the data with negative a in the 
Tables [TV[-|V| (notice that for k = 2.25 there exist negative Gibbs free energy solutions with 
such property). 

Like for the type B configurations, we also find the flux tube solutions with one node, 
relevant data is displayed in the right half of Tables |TV|-|V|. Notice that for k — 1 all the flux 
tube solutions with one node have positive Gibbs free energy. However, with increasing k this 
changes. The important thing is the position of the thermodynamic critical value h c relative 
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to the "limiting points" for h ext of the flux tube families, at which the magnetization and 
the Gibbs free energy vanishes. We have determined numerically the values of such limiting 
points which depend on the value of k and are located, on the h ex t 

axis, at A, ^A, i^A, yA, . . .. 

In the next section we explain the values of these limiting points with the help of linearization 
analysis. 

The reader should compare the magnetization curve and the (AQ/h ext ) energy/field plots 
for k = 1.0, 2.25, Figures [|-|5|, which demonstrate that the "motion" of the limiting points to 
the right of h c is accompanied by creation of flux tubes with one node which have negative 
Gibbs free energy. 

Both, the type B solutions and the flux tube solutions, have a well defined Gibbs free 
energy for an infinite sample. It is natural to compare them. Of course, this must be done 
in a correct way: one should compare energies of configurations with the same values of 
external magnetic field h ext . Using our data, we can display the Gibbs free energy AQ as a 
function of h cxt . These functions for the type B and the flux tubes are given on the Figure 
[5], showing that the flux tube configurations are energetically more preferable. 

6. LINEARIZED SYSTEM AND CRITICAL MAGNETIC FIELDS 

The best way to understand the structure of type B and flux tube solutions in the limit 
of vanishing magnetization M — > and the Gibbs free energy AQ — > is to study the 
linearized Ginzburg-Landau equations. 

Let us consider, in the spirit of |]3|j5|-[T0l , the system ( |2.7|) -(|2~8|) in the situation when 
the square of the scalar field tp 2 is negligibly small. Mathematically this means that, in the 
lowest order, one drops out the terms containing (p 2 in ( |2.7| ) and ( |2.8| ). We then immediately 
notice that such a linearized system 

r 2 f" - r f = 0, (6.1) 

rV' + rV = <p(f 2 - Ar 2 ), (6.2) 

is a consequence of any of the first order systems 
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-/' - sX = 0, (6.3) 
r 

if/ + s-fif = 0, (6.4) 
r 

where s = ±1 is the sign factor. 

This system is straightforwardly integrated and yields 

f(r) = N + s^r\ (6.5) 

tf>(r) = (p r~ sN exp ^— > (6-6) 

where N and y?o are integration constants. In the linear approximation, the constant N is 
arbitrary, but as we can see from the analysis of the complete self- consistent system at the 
origin (|2.18|) -( p.21|) , this constant should be either or ±1, ±2, .... In order to have a regular 



behavior of (|6.6j ) at r = 0, one should choose the sign of N in such a way that sN < 0. 
Notice that the potential (|6.5| ) describes an homogeneous constant magnetic field h = sX, 
and thus the value of s shows its direction (up or down along the z axis). 

It is easy to check that all the solutions ( |6.5| ) and ( |6.6|) , for arbitrary values of integration 



constants, have the same (zero) energy integral computed for the linearized system Q6.1 )- 
( p.2| ). For N = the field ( |6.6| ) evidently describes the linearized type B solution, while 
for iV = 1 this is a linearized flux tube solution. As we see, the linearized solutions are 
energetically equivalent. However, the numerical results (see Figure |5|) definitely show that 
self-consistent flux tubes are energetically more preferable than the type B configurations. 

The first correction to the magnetic field is easily computed. One must now take the 
complete system, and consider the first equation fl2.7|) in the form (|2.11|) where the right 
hand side is constructed from the lowest order configurations (|6.5|) and (|6.6|) . Since these 
satisfy (|6^ ), we find from it \ jip 2 = — |s((^ 2 )', and hence ( |2.11| ) is immediately integrated, 
yielding for the magnetic field 

h = a(\- . (6.7) 

As we see, an important role is played here by the normalization of the linearized solution, 
i.e. by the constant (p . 
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Further insight can be obtained also directly from the analysis of the second order lin- 
earized system (|6.1|) - (|OD . Indeed, integration of flS.ip is straightforward, giving 

f = N+K\ (6.8) 

where N and h are integration constants with the latter representing the value of an homo- 
geneous constant magnetic field. After substituting (|6.8j ) into ( |6.2|) we find a Schrodinger 
type of equation for if with the potential of a circular oscillator. Regular solutions exist only 
when 

sX 

h= l + 2n + siV+|iV|' (6 ' 9) 

where n = 0,1,2, ... . Corresponding eigenfunctions ip Ui N are given in terms of the Laguerre 
polynomials, with n equal to the number of zeros (nodes). Let us introduce the notation 

/l * := 2ifeTT' fc = ' 1 ' 2 '---- ( 6 - 10 ) 



It is easy to see that the maximal eigenvalue ( |6.9| ), h = sX = sho, is achieved for n = and 
sN = —\N\, and the scalar field is then described exactly by ( |6.6| ). This maximal eigenvalue 
is precisely the second critical field h C2 = X = k 2 . The rest of eigenvalues also have clear 
physical meaning: these define the values of the external magnetic field at which the exact 
type B and flux tube solutions become "linearizable" and thus disappear. Looking at the 
Tables |TV|-|V], Pl-jOl, we find the complete agreement with the above linearization analysis. 



Indeed, the flux tube (without nodes) and the type B (without nodes) configurations have 
the limit magnetic field values hi = A/3 and h Q = h C2 , while the flux tubes and type B 
solutions with one-node "live" between h,2 = A/5 and hi = A/3. In general, the family of 
solutions with k nodes have the limiting points hk+i and h^. If h c belongs to the interval 
[hk+i, hk], then hk+i < h ext < h%. for all solutions in this family. However when h c does not 
belong to this interval, then [hk+i, hk] is extended up to h c . 

This linearization analysis clearly supports the existence of the flux tube type solutions. 
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7. OSCILLATING SOLUTIONS 



Let us now consider weaker conditions at infinity: potential still satisfy (3J.), however 
instead of ( |3.2|) we require regularity and flniteness of the scalar field. From the physical 
interpretation of <p as the "density of superconducting electrons" one concludes that \p>\ < 1 
for all values of the radial coordinate r. 

Qualitatively, one can understand the behavior of a scalar field for large r as follows. 
When in (|2.8| ) the potential f 2 and the scalar field p 2 nonlinear terms become small enough, 
one is left with the linearized equation 

<p" + V + \<p = 0, (7.1) 

r 

which has the Bessel function as a solution (p = J (y/\r) = Jq[kt). Such an asymptotic 
behavior is confirmed by direct numerical integration, see Figure |7| and Table [VI|. 

To the best of our knowledge, this type of solution was never reported in the literature. 
It is interesting to find out, what physics corresponds to it. One interpretation is that these 
new oscillating solutions are unstable configurations preceding to the completely formed 
Abrikosov vortex state, they appear when the external magnetic field is switched on and 
reaches H Cl . We may draw attention to the following remarks in an experimental research 
paper: "Consistent values of the magnetization were obtained for fields just above H cl only 
after the sample had been moved between the coils a number of times. The change in the 
magnetization of a sample upon a slight increase (or decrease) of the field is very dependent 
on the fact that the sample has been jarred as it is pulled between the measuring coils, and 
a final, steady-state value of the magnetization is sometimes obtained only after 10 or 20 
sample translations. It is as if vibration assists the flux movement into or out of the sample. 
Hence, all of the data reported below refer to the final steady state of the magnetization; 



that is, further sample motion would produce no further change" [24 



For the type A conditions at the origin, the (selected) results of numerical integration are 
given in Table |VT[ In general, oscillating solutions exist only for initial values (a, b) below 
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the ones (a*,b*) of the vortex solution (first line of Table |V^) . The flux for the displayed 
solutions is always Fj (2tt) = — 1. The values here are given for k = 1.0. From Figure [7|, one 
recognizes that for b close to b* the scalar field reaches almost the value 1, i.e. a complete 
superconducting state. 

Due to a not so quick decay of the scalar field at infinity, approximately if ~ cos(Kr)/y / r, 
the energy of an oscillating solution is infinite for an infinite sample, but for a real finite 
cylindrical sample it is finite, although larger than the energy of a vortex configuration. 
Notice, that an oscillating character of such solutions may resemble an "intermediate" su- 
perconducting state with coexisting normal and superconducting regions f6j. However a 
considerable difference is that the magnetic field penetrates only at the center, exactly like 
in a vortex case. Moreover, the magnetic flux is quantized in a precisely same manner as 
for vortices, which is immediately seen after using ( |2.18| ) and f{oo) = in (|2.5|) , flux is 
Fj (27r) = — N. Notice also that the magnetic energy of oscillating solutions is always finite 
even for an infinite sample. 

Oscillating solutions exist also for the type B conditions. However, the numerical analysis 
revealed that in this case the parameter a must vanish and hence the magnetic field is 
completely absent. Nevertheless, the scalar field configuration is nontrivial. In fact, one is 
left then with the nonlinear scalar field equation 

<f" + -<p' + \(cp — <^ 3 ) = 0, (7.2) 



which in the limit of r — > oo, when ip approaches 0, reduces to the linearized equation ( [7.1|) . 



It is evident, that it is enough to find explicitly only a solution <fi(r) for the case A = k = 
1. For an arbitrary value of the coupling constant the solution is then given by ip = ifi(^/\r). 
In particular, this defines positions of extrema and zeros of <p(r) for all values of A from that 
of (fi(r). The latter evidently depends only on the value of the parameter b in ( |2.21[ ) . 
Numerical results are displayed in Figure [7[ 

As we already mentioned, for an infinite sample both energy and the Gibbs free energy 
are divergent for oscillating solutions. However, for a finite cylinder we discover convergent 



results. As it is well known, the boundary conditions in the Ginzburg-Landau theory require 
vanishing of derivative (p'\ r=R = on a cylinder's surface r = R. Using ( |7.2| ), one then finds 
from ( |2.16| ) that for such zero magnetic field oscillating solutions the Gibbs free energy is 



always negative, 



AQ = [ R drr<p 4 (r). (7.3) 

2 Jo 



A curious conclusion is thus that for a finite sample in absence of an external magnetic field 
an oscillating state is energetically more preferable than a purely normal state. 

We have calculated the Gibbs free energy for the oscillating solutions in finite samples. 
A boundary can be placed at any of the positions of extrema of the solutions, and the 
numerical results are displayed in the lower part of Table [VI|. We present explicitly only the 
case A = 1, while for an arbitrary A, the relevant data are easily obtained from Table [V| by 
replacing (R, AQ) with (R/V\,AQ/X) (cf. (p()). 

Oscillating scalar field solutions appear also if one takes instead of an electromagnetic 
theory a general relativistic gravitational theory |]25|| . The solutions for this Einstein-scalar- 
field theory describe a dark halo of galaxies or galaxy clusters, respectively. The oscillating 
behavior of the scalar field can be removed simply by adding a mass term for the scalar 



field potential. In this case one speaks about boson star solutions p6j which have some 
characteristics similar to the neutron stars but also decisive differences ||27|| . These boson 
stars could be formed in the very early universe from Higgs or axion particles. That stable 
configurations of these boson stars can exist was investigated with the help of the catastrophe 
theory JJ8 . 



8. DISCUSSION AND CONCLUSION 



In each class of solutions the decisive role is played by the values of the parameters (a, b) 
which appear in the regularity conditions at the origin fl2.18| )-( |2.19D and ( |2.2U| )-( |2\21| ). All 



the solutions are obtained after a "fine tuning" of these parameters. It is worthwhile to 
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draw a kind of a "phase diagram" on the (a, b) plane which shows explicitly the domains of 
existence for different solutions. Since the vortices, flux tubes and the oscillating solutions 
all belong to the type A regularity conditions ( |2.18| )-( p.l9| ), we can display them on the same 
(a, o)-plane, see Figure |6| for different values of k. The encircled dots denote the "position" of 
a vortex solution while each curve represents a complete family of a flux tube or oscillating 
solutions for a fixed n. Notice that all curves end on the a-axis (b = 0) at the points which 

correspond to the half of the relevant limit magnetic field values, i.e., a = k — 0, 1, 

Curves which represent flux tubes with increasing number of nodes are concentrating in the 
close neighborhood of the "oscillating" curve which seem to indicate that the oscillating 
solutions are unstable and a small perturbation may cause their decay into a nearby flux 
tube with a finite number of nodes. When moving along any flux tube curve away from the 
a-axis, one inevitably hits the vortex dot, where magnetization diverges. 

Technically, it is impossible (because of the limitations on numerical precision) to make 
integration for the parameters (a, 6) in the close vicinity of a vortex. Thus, from the data 
which we obtained, it is not clear what is the limiting value of an external magnetic field 
to which all the flux tube configurations approach when (a, 6) are coming closer and closer 
to the vortex parameters (a*, 6*). We can see however (cf. Tables |TV|-|V|), that such a limit 
is close to the thermodynamic critical field h c for each k. The following simple argument 
demonstrates that in fact such a limit is equal to h c . Let us formally compare the values of 
the Gibbs free energy for a vortex and for a flux tube solution. We find that these are equal 
when 

S Y - S FT = 2irH cxt (l ~ J™ drr h FT ) , (8.1) 

where the subscripts v and ft denote the vortex and flux tube variables, respectively. The 
right and the left hand sides are both formally divergent, but comparing the leading terms, 



one can use (|8.1|) and (p.6|) to find [noticing that after a certain finite value of r one has 



h F T = h ext and ip FT = 0] h ext = (h 2 cyit + f ) / (2/i cxt ), from which 
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h ext = ]J- = h c . (8.2) 

Below (above) this value, the pure vortices (the flux tubes) are energetically more preferable. 

Summarizing, in this paper we present the numerical solutions of the cylindrically sym- 
metric Ginzburg-Landau equations. Besides the well known vortex configurations with finite 
energy we find new solutions (we call them type B and the flux tube solutions) which have 
finite Gibbs free energy. Direct numerical integration reveals many interesting properties of 
these solutions. One of the most important points is perhaps the clarification of the meaning 
and value of the upper critical field h C2 . Contrary to what is usually claimed in the literature 
(with an exception of [[|), h C2 by no means denotes the magnetic field below which the vortex 
becomes more energetically preferable than the normal state. Instead, as we demonstrated, 
h C2 = A is the value of an external magnetic field at which the type B solutions and the 
flux tubes have zero Gibbs free energy. Below it, for h cxt < h C2 , AQ is negative for both 
flux tubes and type B configurations. The analysis of linearized Ginzburg-Landau equations 
near h C2 which usually (and incorrectly) is described in the literature (see, e.g., [[|0-f§|lO| ) 
as relevant to vortices, in fact is the linearization of flux tubes and type B solutions. Our 
results show that the flux tube solutions without node remain the most energetically prefer- 
able from h C2 down to the thermodynamic critical field h c , after which the vortices become 
energetically more preferable and such a vortex state ends at the lower critical field h ci . We 
find it convenient to depict this observation on Figure 

It is worthwhile to stress, that our results do not contradict the previous knowledge 
about the mixed state in the type II superconductors. On the contrary, they again support 
the significance of such a fundamental structure as a vortex: notice that, after all, one can 
interpret a flux tube solution as a vortex "surrounded" by a type B configuration. However, 
in our opinion, the flux tubes provide us with a new understanding that the mixed state 
reveals a rich structure in which a "pure vortex" is only part of the whole picture valid near 
h ci . We are convinced that a correct transition from such a pure vortex state to the normal 
state (starting at h c up to h C2 ) can only be correctly described with the help of the flux 
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tube and the type B solutions. As is well known, the Ginzburg-Landau theory is only an 
approximation (valid near the critical temperature of superconducting phase transition) to 
the underlying microscopic Gorkov theory within the BCS scheme. The isolated vortices and 



the vortex lattice structures are discussed in the broader aspects in the recent review |29| . It 
is worthwhile to mention the success in theoretical construction of a vortex lattice solution 
for the Gorkov equations [30[|. The work is now in progress aiming at a generalizing our 



isolated flux tube solutions to the lattice type structures. Physically, it would be interesting 
to study a possibility of relating the new solutions to the problem of the origin of the so called 
irreversible line on the phase diagram for the high-temperature superconductors (currently 
discussed explanation of which is the phenomenon of the vortex lattice "melting" |pT|). 

APPENDIX: CORRESPONDENCE WITH THE GINZBURG-LANDAU 

NOTATION 

In the Ginzburg-Landau theory of superconductivity [jT| , the scalar field p is interpreted 
as the "order parameter" with the square describing the "density of superconducting elec- 
trons", n s = \ip\ 2 . The potential is usually written in the form 

V GL = o#| 2 + ^| 4 , (A.l) 

with constant parameters a < 0,(3 > 0. Their physical meaning is clarified by the follow- 
ing quantities they define: the thermodynamic critical magnetic field for a bulk supercon- 
ductor H 2 := a 2 /j3; the equilibrium density of superconducting electrons \<Poo\ '■= \ce\/P; 



the order parameter coherence length £ := l/w2|a|; the magnetic field penetration length 



5 := J(3/\a\. (We are using the units in which the mass and the charge of the electron is 
equal one). Of particular importance is the ratio of two lengths 

k:= 6 -. (A.2) 



Comparing ( |2.2j ) and ( |A.1| ), we find the relation between our and the Ginzburg-Landau 
notations: 
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hence, in our notation, we have 




(A.4) 



(A.5) 



Technically, there are also other notational differences: in the literature on type II super- 
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FIGURES 

FIG. 1. A typical vortex solution. At the center, the magnetic field h is maximal. The limit 
of the scalar field (p is 1, which means physically a completely superconducting state. With our 
gauge choice the potential / is vanishing at infinity. 

FIG. 2. Type B solutions without node (above) and with one node (below) for k = 1.0. In both 
cases the scalar field ip has initial values </?(0) = 0.9, 0.6, 0.3 and magnetic potential /(0) = 0.0. 

FIG. 3. Flux tube solutions without node for k = 0.5, 1.0, 1.5 with the same magnetization 
M/(2tt) = —6.0 and positive Gibbs free energy (above), and the magnetization M/{2tt) = —4.0 
and negative Gibbs free energy (below). The scalar field of a flux tube reaches a maximum before 
going to zero at infinity (and not to 1, as in the case of vortices). The magnetic potential / (not 
shown) starts at 1, is followed by a minimum and then grows with an asymptotic r 2 behavior, 
so that the magnetic field h is asymptotically constant. Note that for k = 0.5 all solutions have 
positive Gibbs free energy, cf. Fig. |I[ 

FIG. 4. Magnetization curves for type B and A solutions for different k. Above, the type B solu- 
tions without node (right curve for each k), and with one node (left curve). For k = 2.25, we notice 
a change of the sign of the Gibbs free energy within a solution family, while for k = 0.5, 1.0, 1.5 we 
find only positive values of Q. Below, for the flux tubes without node for each value k = 0.5, 1.0, 1.5, 
the corresponding limit values h±, hi, h c , and h C2 are drawn. Solid (broken) lines denote negative 
(positive) Gibbs free energy. In all cases the magnetization diverges at h c . For k = 0.5 the point 
h c does not lie within the two limits h\ and h C2 so that one can find also solutions up to h c . 

FIG. 5. AG against the external magnetic field for k = 1.0 and k = 2.25 for the flux tube and 
the type B solutions. In each case the flux tubes have lower values of AQ than the type B solutions. 
Hence, flux tubes are energetically more preferable. 
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FIG. 6. (a, b) diagrams for k = 1.0, 2.25, 5.0. The big dot describes in each case the correspond- 
ing vortex solution. The drawn, broken and dotted lines represent the flux tubes without node, the 
flux tubes with one node, and the oscillating solutions, respectively. The type B solutions cannot 
be compared in these diagrams because they have different initial values. 

FIG. 7. The oscillating solutions for type A (above) and type B (below) initial conditions. 
For an initial value (here, b = 0.777) very near to the one of the vortex, the scalar field closely 
approaches 1 before the oscillation starts. Near the origin, the vortex can be recognized. No type 
A oscillating solution with an initial value above that of the vortex can be found. For a type B 
oscillating solution, the magnetic field vanishes, while the square of the scalar field, the density of 
superconducting electrons, oscillates without any external magnetic field, and every finite sample 
has a negative Gibbs free energy. 

FIG. 8. The general diagram for different solutions. The vortex state is energetically most 
preferable between H Cl and H c , while above H c , the flux tube configurations replace them. Hence, 
H c gains the following physical meaning for a type II superconductor: Above H c , isolated vortices 
come into contact with the external magnetic field and the flux tube solutions are constructed. 
Energetically, the most preferable solution is a flux tube without node. But also the flux tubes 
with nodes do exist there, if k is large enough. Near H C2 , only the flux tube without node exists. Of 
course, the flux tubes can 'live' also in the vortex and the Meissner state but they are energetically 
less preferable. In this way, one finds a rich fine structure of a superconductor's mixed state. 



Figures 1-8 are not enclosed. See instead the Addendum for complete set of figures. 
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TABLES 



TABLE I. Vortex solutions 







N = 


1 






N = 


2 




K 


a 


b 


£M 


V2H C 


a 


b 


£M 


V2H C 


2.00 


-0.46614 


1.35218 


0.38816 


0.38816 


-0.52764 


1.01524 


0.88928 


0.44464 


1.75 


-0.43108 


1.20961 


0.47922 


0.41932 


-0.48000 


0.82862 


1.07872 


0.47194 


1.50 


-0.39344 


1.06659 


0.61094 


0.45821 


-0.42993 


0.65857 


1.34754 


0.50532 


1.25 


-0.35273 


0.92273 


0.81402 


0.50876 


-0.37708 


0.50680 


1.75284 


0.54776 


1.00 


-0.30828 


0.77735 


1.15676 


0.57838 


-0.32096 


0.37219 


2.41915 


0.60478 


0.75 


-0.25904 


0.62916 


1.82190 


0.68321 


-0.26080 


0.25461 


3.67126 


0.68836 


1/V2 


-0.25000 


0.60328 


1.99999 


0.70710 


-0.25000 


0.23614 


3.99999 


0.70710 


0.50 


-0.20314 


0.47525 


3.47164 


0.86791 


-0.19527 


0.15401 


6.64856 


0.83107 



TABLE II. Type B solutions for n = 1.0. 





Without node 




With one 


node 




b 


a 


M/(2vr) 


AG [V] 


hext 


a 


M/(2vr) 


AG [V] 


h e xt 


1 - 


io- 16 


2.62 IO -13 


-319.831 


-3.336 


0.712 


1.88 IO" 14 


-399.137 


22.180 


0.672 


1 - 


io- 12 


2.80 10~ 9 


-150.422 


-2.301 


0.714 


2.62 10~ 10 


-201.853 


14.548 


0.659 


1 - 


10~ 8 


1.79 10~ 6 


-69.240 


-1.582 


0.718 


2.20 10~ 7 


-102.176 


9.272 


0.641 


1 - 


io- 4 


1.11 10~ 3 


-18.774 


-0.846 


0.728 


2.11 10~ 4 


-34.388 


4.108 


0.597 


0.9 




0.1466 


-1.412 


-0.169 


0.800 


0.0439 


-5.299 


0.610 


0.477 


0.7 




0.3030 


-0.430 


-0.042 


0.878 


0.0957 


-2.376 


0.188 


0.420 


0.5 




0.4033 


-0.160 


-0.009 


0.937 


0.1307 


-1.126 


0.051 


0.380 


0.3 




0.4659 


-0.048 


-0.001 


0.977 


0.1536 


-0.401 


0.007 


0.351 


0.1 




0.4962 


-0.005 


-1.25 10~ 5 


0.997 


0.1652 


-0.044 


9.30 10~ 5 


0.335 


0.0 




0.5000 


0.000 


0.000 


1.000 


0.1666 


0.000 


0.000 


0.333 
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TABLE III. Type B solutions for re = 2.25. 





Without node 




With one 


node 






b 


a 


M/(2vr) 


AG [vr/i 2 




h cx t 


a 


M/(2tt) 


a /-i r 

AG [tt/j, 2 






1 - 


io- 8 


3.96 10~ 4 


-66.111 


-2.249 




1.670 


3.96 10~ 4 


-70.821 


0.526 




1.534 


1 - 


io- 4 


0.0396 


-14.189 


-1.107 




1.765 


0.0337 


-17.351 


-4.34 10- 


9 


1.495 


1 - 


10~ 2 


0.3836 


-2.443 


-0.428 




2.069 


0.2154 


-5.055 


-0.101 




1.486 


0.9 




1.1088 


-0.402 


-9.99 10~ 


-2 


2.892 


0.4756 


-1.673 


-3.78 10- 


o 

-z 


1.544 


0.7 




1.7810 


-9.67 10~ 2 


-1.83 1(T 


2 


3.865 


0.6640 


-0.621 


-8.38 10- 


■3 

O 


1.613 


0.5 




2.1691 


-3.33 10~ 2 


-3.52 IO" 


-3 


4.475 


u. / oy t 


— u.zoy 


1 79 i n~ 


-3 




0.3 




2.4044 


-9.81 10~ 3 


-3.90 10- 


4 


4.855 


0.8147 


-8.41 10~ 2 


-1.97 10- 


-4 


1.675 


0.1 




2.5173 


-9.98 10~ 4 


-4.48 10- 


-6 


5.039 


0.8405 


-8.93 10~ 3 


-2.30 10- 


-6 


1.686 


0.0 




2.53125 


0.000 


0.000 




5.0625 


0.84375 


0.000 


0.000 




1.6875 



TABLE IV. Flux-tube solutions for re = 1.00 



Without node 


With one node 


a 


b 


M/(2vr) 


AG M 


h cxt 


a 


b 


M/(2n) 


AG M 


h cxt 


0.16666 


0.00192 


-0.0001 


1.01 10~ 9 


0.3333 


0.09500 


0.05752 


-0.1982 


0.0012 


0.2061 


0.09149 


0.24083 


-0.9999 


0.1087 


0.4408 


0.05000 


0.19570 


-1.6964 


0.0968 


0.2559 


0.00770 


0.38158 


-1.9000 


0.3749 


0.5333 


0.00000 


0.29595 


-3.0598 


0.3264 


0.3041 


-0.00172 


0.39592 


-2.0000 


0.4109 


0.5426 


-0.10000 


0.46591 


-5.5963 


1.0575 


0.3892 


-0.29562 


0.76607 


-12.0000 


1.8288 


0.7304 


-0.28500 


0.74590 


-15.4375 


3.8649 


0.5416 


-0.30818 


0.77728 


-40.0000 


1.3970 


0.7216 


-0.30825 


0.77732 


-56.9794 


8.4137 


0.6195 


-0.30828 


0.77735 


-154.2204 


0.2883 


0.7145 


0.30840 


0.77735 


-39.2231 


4.6106 


0.6053 


0.30828 


0.77735 


-151.5651 


-3.0233 


0.7130 


0.30000 


0.74194 


-10.8059 


1.5583 


0.5326 


0.39637 


0.72724 


-1.5999 


-0.2007 


0.8270 


0.18000 


0.23311 


-0.7157 


0.0176 


0.3582 


0.49998 


0.00999 


-0.0001 


-2.49 10~ 9 


0.9999 


0.16670 


0.01173 


-0.0018 


1.20 10~ 7 


0.3333 
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TABLE V. Flux-tube solutions for re = 2.25 



Without node 


With one node 


a 


b 


M/(2tt) 


AG [TTfi 2 ] 


hext 


a 


b 


M/(2ir) 


AG [V] 


hext 


0.84325 


0.03080 


-0.0010 


0.26 10~ 7 


1.6876 


0.50620 


0.01167 


-0.0003 


0.22 10~ 8 


1.0125 


0.43741 


0.87389 


-1.0000 


0.0146 


1.7754 


0.50000 


0.13056 


-0.0415 


0.0000 


1.0168 


0.00000 


1.23234 


-3.0166 


0.0421 


1.8179 


0.00000 


1.15113 


-3.7504 


0.1824 


1.2955 


-0.01143 


1.23977 


-3.1000 


0.0424 


1.8181 


-0.10000 


1.24623 


-4.8332 


0.2573 


1.3399 


-0.33764 


1.42257 


-8.0000 


-0.0170 


1.7936 


-0.45000 


1.47366 


-17.7012 


0.7550 


1.4838 


-0.49780 


1.49403 


-55.1538 


-1.0370 


1.6803 


-0.48750 


1.48977 


-30.9911 


0.9853 


1.5089 


0.50045 


1.49512 


-49.9903 


-2.3225 


1.6823 


0.50000 


1.49494 


-59.2582 


0.1456 


1.5301 


0.70035 


1.56827 


-5.0000 


-0.8568 


1.9461 


0.55000 


1.51459 


-14.7491 


-0.2253 


1.5043 


1.13578 


1.65419 


-1.2000 


-0.3092 


2.5200 


0.70000 


1.52843 


-4.5843 


-0.0988 


1.5491 


2.53100 


0.03330 


-0.0000 


-4.27 10~ 9 


5.0620 


0.84373 


0.03177 


-0.0005 


0.0000 


1.6874 



TABLE VI. Oscillating solutions for re = 1.0 (R gives the positions of extrema). 



Type A 


a 


b 


R (1st max) 






-0.30828 


0.77735 


oo 


0.24523 




-0.30700 


0.77500 


3.36848 


0.24328 


F/(2n) = -1 


-0.11438 


0.40000 


1.85499 


0.04429 




-0.03500 


0.20000 


1.83312 


0.00579 




-0.00246 


0.05000 


1.84246 


0.00005 


TypeB 


R (min) 


QM 


R (max) 


QW 2 \ 


{a = 0.0, 


3.84 


-2.484 10~ 5 


7.02 


-2.978 10~ 5 


b = 0.1} 


98.2 


-5.009 10~ 5 


208.1 


-5.583 10~ 5 


{a = 0.0, 


5.13 


-0.54496 


8.48 


-0.62000 


b = 0.9} 


295.5 


-1.08209 


198.0 


-1.03113 
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ADDENDUM: COMPLETE SET OF FIGURES AND TABLES 

FIG. 1. A typical vortex solution. At the center, the magnetic field h is maximal. The limit 
of the scalar field ip is 1, which means physically a completely superconducting state. With our 
gauge choice the potential is vanishing at infinity. 

FIG. 2. The scalar field ip and the normalized magnetic field h/(2a) for several values of n. 
From these curves, one can read off the penetration lengths 5 ca i of the scalar field and £ ca z of the 
magnetic field which are determined at the 1/e level. 

FIG. 3. The energy per unit length for vortices with flux F/(2ir) = —1 (N = 1) and 
F/(2ir) = —2 (N = 2). For comparison, the double energy value of the 1-vortex is also given. 
One recognizes clearly that for k > 1/ \/2 (the vertical line) a 1-vortex is energetically more prefer- 
able than a 2-vortex. 

FIG. 4. The comparison of the formal parameter k, which appears in the Ginzburg-Landau 
system, and of the ratio 5 C ai/£,cai calculated from the solution of the Ginzburg-Landau system. 
Near k = 1.0, we find an essential deviation for F/(2ir) = —1, —2. 

FIG. 5. Type B solutions without node for k = 1.0: The scalar field (p with initial values 
(p(0) = 0.9,0.6,0.3 and /(0) = 0.0. 

FIG. 6. The same type B solutions without node for k = 1.0 as in Fig. ||: The potential / and 
the magnetic field h. 

FIG. 7. Type B solutions without node with very high magnetization. The scalar field (p has 
initial values which are almost 1 (see Table V of Addendum). With increasing magnetization the 
scalar field has a wider central core where it is almost 1, i.e. full superconductivity at the center. 
When the scalar field becomes zero, the magnetic field achieves a constant value. 

FIG. 8. Type B solutions with one node for k = 1.0: The scalar field cp with initial values 
<p(Q) = 0.9,0.6,0.3 and /(0) = 0.0. 
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FIG. 9. The same type B solutions with one node for k = 1.0 as in Fig. The potential / and 
the magnetic field h. At zeros of the scalar field, the magnetic field has a step-like behavior. For 
the constant magnetic field, the potential grows as r 2 . 

FIG. 10. Flux tube solutions without node: The scalar field ip for n = 0.5, 1.0, 1.5 with the 
same magnetization M/(2w) = 6.0. The scalar field of a flux tube produces a maximum before 
going to zero (and not to 1, as in the case of vortices). 

FIG. 11. Flux tube solutions without node: The potential / and the magnetic field h for 
k = 0.5, 1.0, 1.5 with the same magnetization M/(2n) = 6.0. The potential starts at 1, is followed 
by a minimum and then grows with an asymptotic r 2 behavior, so that the magnetic field is 
constant. 

FIG. 12. Flux tube solutions without node: Configurations with external magnetic field near 
h c . Both solutions are close to the vortex solution for k = 1.0 which 'lies' between these two flux 
tube solutions; cf. Table XIII of Addendum. The dotted curve has a negative Gibbs free energy, 
while the drawn curve has positive Gibbs free energy. Coming closer and closer to the vortex, the 
magnetization increases without a limit, while the scalar field at the center approaches the vortex 
configuration, and simultaneously the interval grows where it is almost 1. Such a limit flux tube 
consists then of a vortex in the center and a type B solution at higher radial values. The initial 
values (a, b) of this solution are close to that of the vortex; cf. Table I of Addendum. 



FIG. 13. The same flux tube solutions without node as in Fig. 12: potential and magnetic field 
configurations. 

FIG. 14. Flux tube solutions with one node: A solution near h c with negative Gibbs free energy. 
The scalar field has one node while the magnetic field has a 'step' at the node's position. Again 
one can describe this solution as a combination of a vortex at the center and a type B solution 
outwards. 
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FIG. 15. Magnetization curves for type B solutions without (right curve) and with one node 
(left curve) for different k. For k = 2.25, we notice a change of the sign of the Gibbs free energy 
within a solution family. For k = 0.5, 1.0, 1.5 we find only positive values of Q. 

FIG. 16. The magnetization curve for flux tube solutions with k = 0.5, 1.0, 1.5 for a scalar 
field without node. For each k, the corresponding limit values hi, ti2, h c , and h C2 is drawn. For 
instance, the curve for k = 1.0 has a branch with a negative Gibbs free energy in the range 
[h c , /i C2 ] = [0.707,1.0] of the external magnetic field and a branch with positive values in the range 
of [hi, h c ] = [1/3, 0.707]. The magnetization diverges at h c . For k = 0.5 the point h c does not lie 
within the two limits hi and h C2 so that one can find also solutions up to h c . The lower figure 
displays in detail how the branches of solutions with positive and negative AG behave near the h c 
limit (k = 1). 

FIG. 17. Magnetization curve for flux tube solutions with k = 1.0 for a scalar field without 
(limits [hi,h C2 ]) and with one node (limits [/i2,/ii])- In both cases the magnetization diverges at 
the external magnetic field h c which is the limiting point for the vortex. 

FIG. 18. AG against the external magnetic field for k = 1.0 for the flux tube and the type B 
solutions. 

FIG. 19. AG against the external magnetic field for k = 2.25 for the flux tubes and the type B 
solutions. The lower figure gives an enlarged view of the flux tube curve with one node, while the 
broken line denotes type B solutions with one node. In each case the flux tubes have lower values 
of AG than the type B solutions. Hence, flux tubes are energetically more preferable. 

FIG. 20. AG against the external magnetic field for k = 5.0 flux tubes without and with one 
node. Both curves limits into h c . 
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FIG. 21. (a,b) diagrams for k = 1.0,2.25,5.0. The big dot describes in each case the corre- 
sponding vortex solution. The drawn, broken and dotted lines represent the flux tubes without 
node, the flux tubes with one node, and the oscillating solutions, respectively. The type B solutions 
cannot be compared in these diagrams because they have different initial values. 

FIG. 22. The oscillating solutions for type A initial conditions. For an initial value (here, 
b = 0.777) very near to the one of the vortex, the scalar field has a part which is almost 1 before 
the oscillation starts. At the center, the vortex can be recognized. No oscillating solution with an 
initial value above that of the vortex can be found. 

FIG. 23. This figure shows what happens with the magnetic field for an oscillating scalar field. 
At the zeros of the scalar field, the magnetic field has a step-like behavior. As in the case of the 
vortex, the magnetic field quickly vanishes, and the flux of the oscillating solution is quantized. 
Here, we have F/(2n) = — 1. 

FIG. 24. The oscillating solution for type B. Here, a solution with a non-vanishing magnetic 
field does not exist. The square of the scalar field, the density of superconducting electrons, 
oscillates without any external magnetic field, and every finite sample has a negative Gibbs free 
energy. 

FIG. 25. The general diagram for different solutions. The vortex state is energetically most 
preferable between H Cl and H c , while above H c , the flux tube configurations replace them. Hence, 
H c gains the following physical meaning for a type II superconductor: Above H c , isolated vortices 
come into contact with the external magnetic field and the flux tube solutions are constructed. 
Energetically, the most preferable solution is a flux tube without node. But also the flux tubes 
with nodes do exist there, if k is large enough. Near H C2 , only the flux tube without node exists. 
Of course, the flux tubes 'live' also in the vortex and the Meissner state but they are energetically 
less preferable. In this way, we find a rich fine structure of a superconductor's mixed state. 
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Figure 4: 
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Figure 5: 
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Figure 8: 
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Figure 9: 
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Figure 10: 
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Figure 11: 
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Figure 12: 
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Figure 13: 
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Figure 14: 





49 



Figure 15: 
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Figure 16: 
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Figure 17: 
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Figure 18: 
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Figure 20: 
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Figure 21: 
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Figure 22: 
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Figure 23: 
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Figure 24: 
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TABLE I. Vortex solutions for N = 1. 



£[1T[J? 



V2H C 



2.50 
2.25 
2.00 
1.90 
1.80 
1.75 
1.70 
1.60 
1.50 
1.40 
1.30 
1.25 
1.20 
1.10 
1.00 
0.95 
0.90 
0.85 
0.80 
0.75 

0.70710678 

0.70 

0.65 

0.60 

0.55 



-0.52985763 
-0.49897503 
-0.46614139 
-0.45240684 
-0.43829401 
-0.43108902 
-0.42378083 
-0.40884163 
-0.39344716 
-0.37756417 
-0.36115464 
-0.35273894 
-0.34417457 
-0.32657239 
-0.30828665 
-0.29886475 
-0.28924269 
-0.27940823 
-0.26934755 
-0.25904488 
-0.25000000 
-0.24848202 
-0.23763777 
-0.22648713 
-0.21500009 



1.63675335 
1.49452061 
1.35218556 
1.29519757 
1.23816023 
1.20961885 
1.18105958 
1.12387797 
1.06659333 
1.00917811 
0.95159772 
0.92273231 
0.89380826 
0.83575334 
0.77735953 
0.74800633 
0.71852956 
0.68891163 
0.65913181 
0.62916556 
0.60328785 
0.59898348 
0.56855009 
0.53782200 
0.50674553 



0.27263936 
0.32221312 
0.38816478 
0.42091387 
0.45839933 
0.47922576 
0.50162817 
0.55190565 
0.61094895 
0.68104249 
0.76529713 
0.81402661 
0.86804706 
0.99548949 
1.1567616 
1.25418337 
1.36584219 
1.49481540 
1.64508815 
1.82190232 
1.9999998 
2.03227505 
2.28578980 
2.59583844 
2.98163991 



0.34079920 
0.36248976 
0.38816478 
0.39986817 
0.41255939 
0.41932254 
0.42638394 
0.44152452 
0.45821171 
0.47672974 
0.49744313 
0.50876663 
0.52082823 
0.54751921 
0.5783808 
0.59573710 
0.61462898 
0.63529654 
0.65803526 
0.68321337 
0.70710678 
0.71129626 
0.74288168 
0.77875153 
0.81995097 
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0.50 -0.20314008 0.47525314 3.47164166 0.86791041 

0.25 -0.13644748 0.30810000 10.65631899 1.33203987 



TABLE II. Vortex solutions for N = 2. 



K 



£{tth 2 



V2H C 



2.00 
1.75 
1.70 
1.60 
1.50 
1.40 
1.30 
1.25 
1.20 
1.10 
1.00 
0.90 
0.80 
0.75 

0.70710678 
0.70 
0.60 
0.50 



-0.52764322 
-0.48000322 
-0.47019357 
-0.45027413 
-0.42993524 
-0.40915269 
-0.38789871 
-0.37708516 
-0.36614139 
-0.34384344 
-0.32096055 
-0.29743900 
-0.27321196 
-0.26080805 
-0.25000000 
-0.24819367 
-0.22226965 
-0.19527916 



1.01524213 
0.82862511 
0.79239817 
0.72409974 
0.65857003 
0.59580159 
0.53578681 
0.50680969 
0.47851821 
0.42398871 
0.37219117 
0.32311960 
0.27676880 
0.25461243 
0.23614630 
0.23313497 
0.19221625 
0.15401329 



0.88928974 
1.07872319 
1.12484629 
1.22772704 
1.34754217 
1.48856910 
1.65645675 
1.75284680 
1.85917290 
2.10795145 
2.41915478 
2.81759368 
3.34249568 
3.67126654 
3.99999979 
4.05933384 
5.08511208 
6.64856450 



0.44464487 
0.47194139 
0.47805967 
0.49109081 
0.50532831 
0.52099918 
0.53834844 
0.54776462 
0.55775187 
0.57968664 
0.60478869 
0.63395857 
0.66849913 
0.68836247 
0.70710674 
0.71038342 
0.76276681 
0.83107056 
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TABLE III. Type B solutions for k = ^/lj2. 



F/(2tt) 



cxt 



1 - 1(T 16 

1 - io- 12 
1 - io- 8 

1 - 1(T 4 
1 - 10~ 3 

0.9 

0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0.0 



0.81463957 10~ 19 

0.49998905 10~ 12 

0.50000041 10~ 8 

0.0000499975 

0.00049975 

0.04750000 

0.09000000 

0.12750000 

0.16000000 

0.18750000 

0.21000000 

0.22750000 

0.24000000 

0.24750000 

0.25000000 



-506.26448051 
-212.66407036 
-97.79618464 
-26.47560174 
-15.55584217 
-2.24829594 
-1.23667262 
-0.76485295 
-0.48591998 
-0.30398178 
-0.18029540 
-0.09603983 
-0.04114477 
-0.01006986 
-0.00000000 



0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 



0.5 
0.5 
0.5 
0.5 
0.5 
0.5 
0.5 
0.5 
0.5 
0.5 
0.5 
0.5 
0.5 
0.5 
0.5 
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TABLE IV. Type B solutions for k = 0.5. 



M/(2vr) 



AG [tt/i 2 



h, 



cxt 



1 - 1(T 16 

1 - io- 12 

1 - 1(T 8 
1 - 1(T 4 

1 - io- 2 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



0.29397206 10~ 29 

0.22785603 10~ 17 

0.11163709 10" 11 

0.57328821 10~ 6 

0.00042699 

0.01112676 

0.02804562 

0.04647666 

0.06469237 

0.08160630 

0.09643545 

0.10858262 

0.11759100 

0.12313071 

0.12500000 



-842.90290064 
-303.48173375 
-140.17105915 
-38.35953571 
-11.04294492 
-3.50056424 
-2.01439886 
-1.29869893 
-0.85736327 
-0.55534931 
-0.33958193 
-0.18551044 
-0.08099744 
-0.02005903 
-0.00000000 



11.741127 

7.0628467 

4.9082144 

2.5230830 

1.32277094 

0.62184922 

0.36480792 

0.21255415 

0.11638777 

0.05742286 

0.02412093 

0.00781468 

0.00157347 

0.00009957 

0.00000000 



0.35182244 
0.35066627 
0.34938890 
0.34539301 
0.33805432 
0.32363276 
0.31151388 
0.29980096 
0.28851017 
0.27797881 
0.26860687 
0.26079463 
0.25490781 
0.25124415 
0.25000000 
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TABLE V. Type B solutions for k = 1.0 



M/(2vr) 



AG [V 



h, 



cxt 



1 - 10~ 16 

1 - io- 12 
1 - io- 8 

1 - 10~ 4 
1 - 10~ 3 

0.9 

0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0.0 



0.2626738 10~ 12 

2.8028373 IO" 9 

0.00000179 

0.00111758 

0.00561929 

0.14665519 

0.23437248 

0.30301143 

0.35848496 

0.40331923 

0.43886076 

0.46590309 

0.48493325 

0.49624589 

0.50000000 



-319.83194287 
-150.42278724 
-69.24007330 
-18.77412288 
-10.98561767 
-1.41282055 
-0.72891143 
-0.43023181 
-0.26377020 
-0.16055448 
-0.09326563 
-0.04892564 
-0.02074118 
-0.00504513 
-0.00000000 



-3.33629143 
-2.30158296 
-1.58213465 
-0.84610904 
-0.65155679 
-0.16978151 
-0.08386589 
-0.04233326 
-0.02052442 
-0.00913650 
-0.00352478 
-0.00106759 
-0.00020468 
-0.00001257 
-0.00000000 



0.71216412 
0.71451434 
0.71804013 
0.72827673 
0.73506133 
0.80012302 
0.84215100 
0.87880681 
0.91069735 
0.93782855 
0.96013136 
0.97754071 
0.99000792 
0.99750049 
1.00000000 
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TABLE VI. Type B solutions for k = 1.5 



a 



M/(2tt) 



h, 



cxt 



1 - io- 12 

1 - 1(T 8 

1 - 10~ 4 

1 - IO" 2 

0.9 

0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0.0 



0.0000009048 

0.00009492 

0.01076082 

0.12487016 

0.43216305 

0.61870985 

0.75504470 

0.86159257 

0.94608432 

1.01229462 

1.06230783 

1.09734361 

1.11811540 

1.12500000 



-120.08566366 
-55.45237424 
-14.51995526 
-3.50311345 
-0.78071937 
-0.36990133 
-0.20814435 
-0.12380358 
-0.07384290 
-0.04230148 
-0.02198078 
-0.00926097 
-0.00224481 
-0.00000000 



-2.90735164 
-2.02045383 
-1.07824205 
-0.51323746 
-0.15951625 
-0.07057877 
-0.03333402 
-0.01546007 
-0.00667256 
-0.00251856 
-0.00075128 
-0.00014259 
-0.00000871 
-0.00000000 



1.08671648 
1.09897771 
1.13651774 
1.22807544 
1.47054176 
1.65414678 
1.80288132 
1.92600329 
2.02721532 
2.10843492 
2.17077409 
2.21490851 
2.24124437 
2.25000000 
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TABLE VII. Type B solutions for k = 2.25 



a 



M/(2vr) 



AG [tt M 2 



h, 



cxt 



1 - 1(T 8 

1 - 1(T 4 

1 - 10~ 2 

0.9 

0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0.0 



0.00039686187 

0.039632414 

0.38362876 

1.1088587 

1.5013300 

1.7810811 

1.9978592 

2.1691771 

2.3032325 

2.4044317 

2.4753075 

2.5173243 

2.5312500 



-66.11152808 
-14.18949823 
-2.44323885 
-0.40297816 
-0.17747571 
-0.09671914 
-0.05652767 
-0.03335736 
-0.01897809 
-0.00981733 
-0.00412457 
-0.00099821 
0.00000000 



-2.2499749 
-1.1073683 
-0.42859875 
-0.099959198 
-0.040526105 
-0.018347512 
-0.0083010393 
-0.0035262489 
-0.0013169381 
-0.00039005489 
-0.000073694017 
-0.0000044888587 
0.00000000 



1.67050641 
1.76576571 
2.06919235 
2.89279937 
3.44436324 
3.86542221 
4.20301812 
4.47522932 
4.69097359 
4.85521398 
4.97087202 
5.03966430 
5.06250000 
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TABLE VIII. Type B solutions with one node for re = 0.5. 



a 



M/(2tt) 



AG [V 



h, 



cxt 



1 - io- 12 

1 - 10~ 8 

1 - 1(T 4 

1 - IO" 2 

0.9 

0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0.0 



0.20777588 10 

0.16968270 10 

1.98235 10~ 8 

0.00003060 

0.00145092 

0.00470204 

0.00925886 

0.01475161 

0.02078605 

0.02691416 

0.03262722 

0.03736276 

0.04053971 

0.04144444 



-19 
-13 



-407.83814571 
-207.47249709 
-70.92516787 
-27.38806803 
-12.10387413 
-8.24390469 
-6.05828898 
-4.49834233 
-3.26222813 
-2.23078265 
-1.36003756 
-0.65666061 
-0.17538241 
-0.00000000 



76.685566 
50.180428 
23.94263827 
11.06437630 
4.79404931 
2.91014633 
1.82769449 
1.10800369 
0.61837208 
0.29897268 
0.11262889 
0.02623761 
0.00185854 
0.00000000 



0.32415859 
0.31279238 
0.28591112 
0.24968540 
0.20705162 
0.18357733 
0.16452267 
0.14740092 
0.13157619 
0.11700553 
0.10401313 
0.09329940 
0.08598007 
0.08333333 
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TABLE IX. Type B solutions with one node for k = 1.0. 



a 



M/(2vr) 



h, 



cxt 



- 10~ 16 

- io- 12 
-lo- 8 

- 10" 4 

- 10~ 3 

-io- 2 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



0.18816282 10 

0.26227586 10 

0.00000022 

0.00021123 

0.00123249 

0.00740425 

0.04397343 

0.07272152 

0.09578531 

0.11488290 

0.13071774 

0.14359881 

0.15364041 

0.16085504 

0.16520978 

0.16666666 



-13 
-9 



-399.13713984 
-201.85338653 
-102.17627279 
-34.38826735 
-22.55431617 
-12.84295319 
-5.29912852 
-3.41688303 
-2.37610307 
-1.66122221 
-1.12628291 
-0.71443290 
-0.40169570 
-0.17911859 
-0.04493592 
-0.00000000 



22.180605 

14.548888 
9.27215886 

4.1088762 
2.86626778 

1.6778879 
0.61000069 
0.33196266 
0.18841580 
0.10272050 
0.05103931 
0.02171670 
0.00713976 
0.00145678 
0.00009304 
0.00000000 



0.67281714 
0.65939048 
0.64105359 
0.59787403 
0.57561429 
0.54090402 
0.47769745 
0.44532585 
0.42035415 
0.39901535 
0.38039878 
0.36441468 
0.35131626 
0.34150658 
0.33540631 
0.33333333 
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TABLE X. Type B solutions with one node for k = 1.5. 



M/(2tt) 



h, 



cxt 



1 - 10~ 16 

1 - io- 12 

1 - 1(T 8 
1 - 10~ 4 

1 - io- 2 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



0.19551488 10~ 8 

0.52845869 IO" 6 

0.48741567 10~ 4 

0.48556580 IO" 2 

0.05185846 

0.16380734 

0.22458591 

0.26661104 

0.29844532 

0.32324161 

0.34247940 

0.35693478 

0.36703569 

0.37301794 

0.37500000 



-253.65744393 
-137.93602705 
-69.41302070 
-22.89693286 
-8.13129020 
-3.07893180 
-1.88505320 
-1.25602039 
-0.84488327 
-0.55284397 
-0.33962509 
-0.18577883 
-0.08108642 
-0.02006596 
-0.00000000 



6.8992671 
4.5873181 
2.7590523 
1.0335935 
0.31314228 
0.08550524 
0.04289325 
0.02314474 
0.01207324 
0.00574979 
0.00234955 
0.00074493 
0.00014762 
0.00000925 
0.00000000 



1.01900652 
1.00509209 
0.98413629 
0.93620705 
0.87891501 
0.83074850 
0.81115382 
0.79687633 
0.78493617 
0.77469908 
0.76608502 
0.75918366 
0.75412788 
0.75103921 
0.75000000 
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TABLE XL Type B solutions with one node for k = 2.25. 



a 



M/(2vr) 



AG [V 



h, 



ext 



1 - io-« 
1 - io~ 4 
1 - io- 2 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



0.0003961775 

0.033788231 

0.21543634 

0.47568332 

0.58990955 

0.66407250 

0.71839505 

0.75977286 

0.79135512 

0.81479746 

0.83103287 

0.84059205 

0.84375000 



-70.82155158 
-17.35196356 
-5.05538676 
-1.67363013 
-0.96900019 
-0.62129733 
-0.40560746 
-0.25909793 
-0.15613863 
-0.08415334 
-0.03634436 
-0.00893709 
0.00000000 



0.52680199 
-0.043408899 
-0.10152137 
-0.037872528 
-0.017383337 
-0.0083894957 
-0.0039511917 
-0.0017251285 
-0.0006568765 
-0.0001972193 
-0.0000376004 
-0.23022948 10~ 5 
0.00000000 



1.53403702 
1.49548933 
1.48666351 
1.54467090 
1.58483248 
1.61373320 
1.63573439 
1.65275343 
1.66581546 
1.67552249 
1.68224194 
1.68619486 
1.68750000 
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TABLE XII. Flux-tube solutions for k = 0.50 



M/(2vr) 



cxt 



0.04153310 

0.03881390 

0.03377073 

0.03187151 

0.02547804 

0.02310788 

0.00943295 

-0.00711859 

-0.02593361 

-0.06611321 

-0.12500000 

-0.14000000 

-0.17000000 

-0.18000000 

-0.19000000 

0.20029645 

0.19100000 

0.17800000 

0.17500000 

0.16000000 

0.15000000 

0.14000000 

0.13250000 

0.12750000 

0.12500000 



0.00482634 
0.02291673 
0.03987468 
0.04510360 
0.06080732 
0.06616722 
0.09464338 
0.12627010 
0.16051533 
0.23111422 
0.33318466 
0.35942324 
0.41279518 
0.43102350 
0.44962804 
0.46875000 
0.44513701 
0.40691601 
0.39713515 
0.34108654 
0.29422019 
0.23346883 
0.16870544 
0.09906230 
0.00100000 



-0.01000000 

-0.20000000 

-0.49999998 

-0.59999997 

-0.89999996 

-1.00000005 

-1.49999992 

-2.00000011 

-2.50000020 

-3.50000033 

-5.25759258 

-5.89146805 

-7.87229590 

-9.05838675 

-11.08773185 

-14.60379550 

-7.67322461 

-4.65587992 

-4.20761305 

-2.52520702 

-1.69236974 

-0.97543913 

-0.48035611 

-0.15975047 

-0.00001600 



0.00001160 
0.00485667 
0.03251648 
0.04786481 
0.11465226 
0.14434866 
0.35425328 
0.67207379 
1.09285215 
2.14659123 
4.0155314 
4.5611306 
5.7877739 
6.2658005 
6.8232515 
2.2039074 
1.4931287 
0.94553126 
0.84563339 
0.43874152 
0.23906195 
0.09420080 
0.02572468 
0.00307150 
0.31997933 10~ 10 



0.08362321 
0.08933301 
0.09914171 
0.10262833 
0.11372799 
0.11763453 
0.13855274 
0.16122362 
0.18459837 
0.22868403 
0.28313110 
0.29539602 
0.31812023 
0.32521571 
0.33221084 
0.33977157 
0.33119462 
0.31967926 
0.31685806 
0.30136833 
0.28941054 
0.27567022 
0.26378724 
0.25485285 
0.25000000 
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TABLE XIII. Flux-tube solutions for re = 1.00 



M/(2vr) 



cxt 



0.16666065 

0.16048691 

0.09149278 

0.00770227 

-0.00172647 

-0.23437054 

-0.28685815 

-0.29562992 

-0.30000000 

-0.30326892 

-0.30818466 

-0.30827686 

-0.30828641 

-0.3082866485 

0.3082866575 

0.308286675 

0.30828685 

0.30828800 

0.30840000 

0.32000000 

0.39637726 

0.44486258 

0.48822938 

0.49875776 

0.49998750 



0.00192455 

0.06238119 

0.24083303 

0.38158876 

0.39592250 

0.70095040 

0.75749117 

0.76607871 

0.77018431 

0.77314720 

0.77728463 

0.77735251 

0.77735936 

0.77735952553 

0.7773595318839 

0.77735954 

0.77735967 

0.77736047 

0.77743589 

0.78196826 

0.72724186 

0.59258754 

0.29977039 

0.09944583 

0.00999947 



-0.00010000 
-0.10000000 
-0.99999998 
-1.90000005 
-2.00000003 
-5.99999872 
-10.00000276 
-12.00000628 
-13.75256885 
-15.99995420 
-40.00000208 
-60.00005351 
-99.99974294 
-154.22042007 
-151.56514834 
-129.08124738 
-100.37162393 
-78.04813009 
-36.78370133 
-8.71898303 
-1.59999998 
-0.60000000 
-0.10000000 
-0.01000000 
-0.00010000 



0.10185308 10~ 8 
0.00103043 
0.10874652 
0.37497658 
0.41094947 
1.54794907 
1.81154073 
1.82881244 
1.82095047 
1.79510039 
1.39704713 
1.12951015 
0.70930334 
0.28830787 
-3.02337416 
-2.41023205 
-2.17454851 
-1.96478896 
-1.46832294 
-0.81916267 
-0.20076202 
-0.05241147 
-0.00225057 
-0.00002473 
-0.000000002499 



0.33334352 
0.34360827 
0.44085382 
0.53338571 
0.54261067 
0.71152198 
0.72966737 
0.73048079 
0.73018163 
0.72931654 
0.72168651 
0.71895537 
0.71624176 
0.71452845 
0.71300061 
0.71510625 
0.71614669 
0.71733231 
0.72196383 
0.74119260 
0.82705431 
0.89933482 
0.97688934 
0.99752029 
0.99997500 
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TABLE XIV. Flux-tube solutions for k = 1.50 



M/(2vr) 



AG [V] 



cxt 



0.37338330 

0.35865108 

0.29017559 

0.20216876 

-0.08153028 

-0.32032191 

-0.37526158 

-0.38396193 

-0.38882706 

-0.39033430 

0.39368037 

0.39942951 

0.40420067 

0.42313776 

0.43745867 

0.46163405 

0.50659366 

0.60827205 

0.65684841 

0.72822112 

0.84790088 

1.04053433 

1.11531973 

1.12401729 

1.12490158 



0.04336260 
0.13882713 
0.32503549 
0.47718137 
0.81864927 
1.02063593 
1.05617319 
1.06126211 
1.06402751 
1.06487144 
1.06672106 
1.06981657 
1.07231309 
1.08161666 
1.08803048 
1.09765879 
1.11151317 
1.12205983 
1.11579485 
1.09126036 
1.00003242 
0.62739169 
0.22190830 
0.07105192 
0.02249687 



-0.01000000 
-0.10000000 
-0.50000000 
-1.00000001 
-2.99999992 
-8.09999184 
-15.00016168 
-19.00079350 
-24.00200915 
-26.99363351 
-43.31969516 
-20.00000303 
-16.00005584 
-10.00000397 
-8.00000030 
-6.00000061 
-4.00000004 
-2.00000014 
-1.50000002 
-1.00000002 
-0.50000000 
-0.10000000 
-0.01000000 
-0.00100000 
-0.00010000 



0.00000786 
0.00077989 
0.01852335 
0.06744250 
0.34444650 
0.58261046 
0.46855172 
0.36778303 
0.24189747 
0.16315885 
-5.99946594 
-1.63731205 
-1.50098647 
-1.22846970 
-1.10316720 
-0.94558871 
-0.73520293 
-0.42922215 
-0.32846827 
-0.21488613 
-0.09110741 
-0.00695891 
-0.00008535 
-0.00000087 
-0.00000001 



0.75176971 
0.76758881 
0.83476960 
0.90884380 
1.07986153 
1.14350498 
1.13288052 
1.12618616 
1.11955789 
1.11607308 
1.00392926 
1.12470204 
1.13328380 
1.15769607 
1.17348368 
1.19913931 
1.24765087 
1.36989865 
1.43530562 
1.53932635 
1.73274206 
2.08404295 
2.23067648 
2.24803495 
2.24980316 
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TABLE XV. Flux-tube solutions for re = 2.25 



M/(2vr) 



h, 



cxt 



0.84325497 

0.79515994 

0.43741045 

0.00233172 

0.00000000 

-0.01143546 

-0.11487304 

-0.33764968 

-0.38602425 

-0.49726408 

-0.49780000 

0.50045166 

0.52081765 

0.58020017 

0.70035444 

1.00000000 

1.13578728 

1.40000000 

1.70971180 

2.12946840 

2.43883772 

2.45577270 

2.50000000 

2.52100000 

2.53100000 



0.03080636 
0.30501597 
0.87389870 
1.23081820 
1.23234317 
1.23977546 
1.30363200 
1.42257539 
1.44534468 
1.49381678 
1.49403737 
1.49512705 
1.50339534 
1.52647271 
1.56827220 
1.63985366 
1.65419047 
1.64269550 
1.54816802 
1.21592253 
0.62758267 
0.56931287 
0.36990731 
0.21281951 
0.03330876 



-0.00100000 
-0.10000000 
-1.00000000 
-3.00000000 
-3.01667010 
-3.10000000 
-4.00000000 
-8.00000000 
-10.00000001 
-50.03000994 
-55.15382089 
-49.99032149 
-20.00000037 
-10.00000000 
-5.00000000 
-1.74625522 
-1.20000000 
-0.62208273 
-0.30000000 
-0.09568099 
-0.01711513 
-0.01380428 
-0.00553393 
-0.00178806 
-0.00004330 



0.26603798 10~ 7 
0.25079340 10~ 3 
0.014681371 
0.042078389 
0.042145730 
0.04244274 
0.04174682 
-0.01708379 
-0.06326788 
-0.97225009 
-1.0370393 
-2.3225981 
-1.6410071 
-1.2232806 
-0.85686163 
-0.42190989 
-0.30921646 
-0.16450278 
-0.07041132 
-0.01331646 
-0.00060838 
-0.00040280 
-6.7727614 10~ 5 
-7.2206406 10~ 6 
-4.2770818 10~ 9 



1.68763471 
1.70038988 
1.77541273 
1.81785677 
1.81791344 
1.81815942 
1.81776491 
1.79367785 
1.78066590 
1.68350270 
1.68033109 
1.68233836 
1.73773679 
1.81341932 
1.94614089 
2.32313287 
2.52009576 
2.94138042 
3.48319881 
4.27157414 
4.87826672 
4.91193731 
5.00006609 
5.04200706 
5.06200000 
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TABLE XVI. Flux-tube solutions for k = 5.00 



M/(2vr) 



h, 



oxt 



4.16596082 

4.09778819 

3.59756383 

2.46720689 

1.12425543 

0.28823276 

0.00000000 

-0.06903770 

-0.10000000 

-0.21388850 

-0.40000000 

-0.50000000 

0.90370250 

0.95902524 

1.05510050 

1.23355635 

1.65490767 

3.19773261 

4.18000000 

5.00000000 

7.14003165 

8.61436325 

9.94894093 

12.00000000 

12.48466309 



0.04810398 
0.47269111 
1.30756971 
2.08399110 
2.58452071 
2.81635033 
2.88809220 
2.90474500 
2.91214944 
2.93905032 
2.98190822 
3.00438855 
3.09142010 
3.10289605 
3.12257405 
3.15828834 
3.23837892 
3.48019659 
3.58538926 
3.63660608 
3.57240809 
3.31817152 
2.88171879 
1.39965274 
0.25004504 



-0.00010000 
-0.01000000 
-0.10000000 
-0.50000000 
-2.00000000 
-4.99999563 
-7.25261516 
-8.00000145 
-8.37263989 
-9.99944370 
-14.06339412 
-17.68306279 
-25.02259657 
-20.00180865 
-14.99986439 
-10.00002257 
-5.00000000 
-1.00000000 
-0.48951133 
-0.29690981 
-0.10000000 
-0.05000000 
-0.02500000 
-0.00349371 
-0.00010000 



-3.4267426 10~ 9 
-3.3077416 10~ 5 
-0.00250402 
-0.02888276 
-0.13149506 
-0.28008789 
-0.36536655 
-0.39067762 
-0.40286046 
-0.45318461 
-0.56348072 
-0.64910298 
-1.3240437 
-1.2461821 
-1.0959379 
-0.90417624 
-0.62843309 
-0.22660722 
-0.13181569 
-0.08628561 
-0.02872954 
-0.01209595 
-0.00440641 
-0.00013615 
-1.2234457 10~ 7 



8.33247708 

8.24994494 

7.65560143 

6.41076481 

5.19549809 

4.61185532 

4.43486702 

4.39333130 

4.37472398 

4.30597683 

4.18974698 

4.12187923 

4.07299171 

4.11784001 

4.22641672 

4.42267078 

4.91311818 

7.16833774 

8.85545777 

10.34339701 

14.40566028 

17.28374629 

19.91846903 

24.00065683 

24.96932678 
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TABLE XVII. Flux-tube solutions with one node for k = 1.00 



M/(2vr) 



h, 



cxt 



0.10000000 

0.09500000 

0.09000000 

0.05000000 

0.02000000 

0.00000000 

-0.10000000 

-0.20000000 

-0.25000000 

-0.28500000 

-0.30500000 

-0.30800000 

-0.30820000 

-0.30825000 

0.30840000 

0.30864000 

0.308645625 

0.30800000 

0.30000000 

0.25000000 

0.20000000 

0.18000000 

0.17000000 

0.16700000 

0.16670000 



0.00000000 
0.05752439 
0.08208260 
0.19570689 
0.25807049 
0.29595722 
0.46591028 
0.62101985 
0.69538079 
0.74590688 
0.77336546 
0.77706223 
0.77727643 
0.77732618 
0.77735374 
0.77650569 
0.77627525 
0.77141194 
0.74194860 
0.57825777 
0.36647786 
0.23311025 
0.11709493 
0.03709171 
0.01173153 



0.00000000 
-0.19829029 
-0.38793945 
-1.69647232 
-2.53551297 
-3.05989492 
-5.59633932 
-8.76927872 
-11.44975608 
-15.43754890 
-25.16144125 
-40.62038139 
-49.74797375 
-56.97949955 
-39.22317054 
-25.20251927 
-24.24790512 
-17.49003861 
-10.80591442 
-4.54518747 
-1.73955332 
-0.71579493 
-0.18374668 
-0.01855901 
-0.00185797 



0.00000000 
0.00122474 
0.00474348 
0.09688239 
0.22211571 
0.32649272 
1.05752348 
2.16007539 
2.96501895 
3.86498975 
5.31339725 
6.96751752 
7.80093365 
8.41376235 
4.61063069 
3.26378469 
3.16433654 
2.41915630 
1.55832296 
0.51406021 
0.09674361 
0.01761512 
0.00119903 
0.00001235 
0.00000012 



0.20000000 
0.20615685 
0.21215358 
0.25596041 
0.28553094 
0.30418516 
0.38929139 
0.46739533 
0.50763559 
0.54164162 
0.57868048 
0.60456234 
0.61383054 
0.61958525 
0.60533567 
0.58397801 
0.58196662 
0.56390083 
0.53261908 
0.45932507 
0.39071118 
0.35821647 
0.33987571 
0.33399871 
0.33339999 
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TABLE XVIII. Flux-tube-solutions with one node for k = 1.50 



M/(2vr) 



AG [V 



ext 



0.22500000 
0.22000000 
0.21000000 
0.20000000 
0.19000000 
0.18000000 
0.17000000 
0.16000000 
0.15000000 
0.14000000 
0.13000000 
0.12000000 
0.11000000 
0.10000000 
0.09000000 
0.08000000 
0.07000000 
0.06000000 
0.05000000 
0.04000000 
0.03000000 
0.02000000 
0.01000000 



0.00000000 
0.08313281 
0.14475946 
0.18785321 
0.22339176 
0.25454495 
0.28275365 
0.30881502 
0.33322364 
0.35631066 
0.37831017 
0.39939456 
0.41969495 
0.43931360 
0.45833199 
0.47681616 
0.49482045 
0.51239011 
0.52956326 
0.54637228 
0.56284492 
0.57900511 
0.59487361 



0.00000000 
-0.08422212 
-0.24910164 
-0.40973461 
-0.56665414 
-0.72032371 
-0.87115056 
-1.01949636 
-1.16568537 
-1.31001117 
-1.45274190 
-1.59412474 
-1.73438975 
-1.87375249 
-2.01241726 
-2.15057901 
-2.28842532 
-2.42613867 
-2.56389763 
-2.70187844 
-2.84025694 
-2.97920954 
-3.11891442 



0.00000000 
0.00019629 
0.00172214 
0.00466970 
0.00894568 
0.01446994 
0.02117316 
0.02899498 
0.03788266 
0.04778985 
0.05867577 
0.07050440 
0.08324393 
0.09686622 
0.11134644 
0.12666273 
0.14279588 
0.15972916 
0.17744808 
0.19594026 
0.21519532 
0.23520479 
0.25596195 



0.45000000 
0.45524174 
0.46553805 
0.47560319 
0.48545641 
0.49511432 
0.50459133 
0.51390008 
0.52305168 
0.53205603 
0.54092196 
0.54965737 
0.55826941 
0.56676455 
0.57514868 
0.58342721 
0.59160507 
0.59968684 
0.60767672 
0.61557866 
0.62339630 
0.63113309 
0.63879225 
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TABLE XIX. Flux-tube solutions with one node for k = 2.25 



M/(2vr) 



AG [vr/x 2 



cxt 



0.50620000 

0.50500000 

0.50000000 

0.10000000 

0.00000000 

-0.10000000 

-0.30000000 

-0.45000000 

-0.48500000 

-0.48750000 

0.50000000 

0.50750000 

0.52000000 

0.55000000 

0.57500000 

0.60000000 

0.65000000 

0.70000000 

0.75000000 

0.80000000 

0.82500000 

0.83500000 

0.84000000 

0.84370000 

0.84373500 



0.01167690 
0.05838531 
0.13056159 
1.03984381 
1.15113625 
1.24623535 
1.39490448 
1.47366214 
1.48872942 
1.48977303 
1.49494165 
1.49800748 
1.50304576 
1.51459409 
1.52330439 
1.53070812 
1.53870849 
1.52843226 
1.47278637 
1.27208198 
0.97248090 
0.71480687 
0.48692825 
0.05799446 
0.03177454 



-0.00033251 
-0.00831132 
-0.04152946 
-2.85901191 
-3.75047509 
-4.83323620 
-8.40318054 
-17.70122111 
-28.96559305 
-30.99115126 
-59.25824667 
-32.01239549 
-22.66014509 
-14.74916477 
-11.58013568 
-9.45207563 
-6.58697807 
-4.58431344 
-2.95284051 
-1.43082804 
-0.63822309 
-0.30403051 
-0.13178990 
-0.00177295 
-0.00053194 



0.22768052 10" 
0.00000142 
0.00003540 
0.12071005 
0.18242766 
0.25732584 
0.46201774 
0.75503469 
0.95568218 
0.98537998 
0.14568520 
-0.12085067 
-0.19695500 
-0.22530268 
-0.21460986 
-0.19508325 
-0.14762030 
-0.09886668 
-0.05379746 
-0.01658306 
-0.00388757 
-0.00095361 
-0.00018702 
-3.50067 10~ 8 
0.00000000 



1.01253467 
1.01336606 
1.01681860 
1.24799509 
1.29551424 
1.33996121 
1.42158269 
1.48389431 
1.50645779 
1.50896690 
1.53019652 
1.51498887 
1.50793798 
1.50431979 
1.50645849 
1.51122284 
1.52662205 
1.54911417 
1.57988712 
1.62390581 
1.65544607 
1.67131472 
1.68025546 
1.68740005 
1.68747000 
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TABLE XX. Flux-tube solutions with one node for k = 5.00 



M/(2tt) 



AG [V 



cxt 



4.16660000 

4.00000000 

3.50000000 

2.90000000 

2.30000000 

1.30000000 

1.20000000 

1.10000000 

1.00000000 

0.90000000 

-0.69000000 

-0.60000000 

-0.50000000 

0.00000000 

0.50000000 

1.00000000 

1.50000000 

2.00000000 

2.25000000 

2.35000000 

2.45000000 

2.47500000 

2.49500000 

2.49750000 

2.49975000 



0.03194924 
1.54205972 
2.74740880 
3.23595155 
3.31897602 
3.17127693 
3.15164667 
3.13166024 
3.11132720 
3.09064825 
3.04608002 
3.02649613 
3.00438828 
2.88783642 
2.75540523 
2.57972937 
2.29453336 
1.76681579 
1.30125018 
1.02401580 
0.60047219 
0.42623100 
0.19120008 
0.13525045 
0.04278467 



-0.00002206 
-0.05963664 
-0.30980879 
-0.85308464 
-1.97604538 
-9.15487152 
-11.20374194 
-14.07329927 
-18.46102014 
-26.43953266 
-37.39455337 
-25.05396233 
-18.71862741 
-7.65473040 
-3.98870419 
-2.15981734 
-1.10535940 
-0.44288754 
-0.20067548 
-0.11597628 
-0.03727405 
-0.01847061 
-0.00366787 
-0.00183231 
-0.00018309 



-1.1760422 10~ 10 
-0.00076581 
-0.01394038 
-0.05940588 
-0.15417819 
-0.52003242 
-0.58424333 
-0.65910757 
-0.75012231 
-0.87154430 
-0.44849527 
-0.36250996 
-0.30676287 
-0.16951423 
-0.09554523 
-0.04797991 
-0.01910902 
-0.00431046 
-0.00102708 
-0.00036295 
-0.00003960 
-0.00000986 
-0.00000039 
-9.8159338 10~ 8 
-9.8120083 10~ 10 



8.33320000 
8.00283544 
7.05068132 
6.00933253 
5.11558078 
4.03145976 
3.95212203 
3.87747347 
3.80669855 
3.73762810 
3.68172157 
3.71708280 
3.74928043 
3.89400437 
4.06148335 
4.26240268 
4.49078655 
4.73836093 
4.86761408 
4.92020750 
4.97328463 
4.98662776 
4.99732323 
4.99866147 
4.99986614 
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TABLE XXI. Oscillating solutions for k 


= 1.0 and flux F/(2?r) = -1. 




a 


b 


£mag [^A* ] 


R (1st max) 


-0.30828665 


0.77735953 


0.24523135 


oo 


-0.30809193504244 


0.77700000 


0.24493325 


4.02791712 


-0.30700707 


0.77500000 


0.24328646 


3.36848331 


-0.26629161 


0.70000000 


0.18786525 


2.21217601 


-0.21314184 


0.60000000 


0.12758654 


1.99170090 


-0.16217784 


0.50000000 


0.08008222 


1.89958297 


-0.11438161 


0.40000000 


0.04429725 


1.85499121 


-0.07120571 


0.30000000 


0.01975218 


1.83620502 


-0.03500354 


0.20000000 


0.00579848 


1.83312370 


-0.0095547963444034 


0.10000000 


0.00058777 


1.83776768 


-0.0024645497256301 


0.05000000 


0.00005105 


1.84246831 



TABLE XXII. Type B oscillating solutions for n = 1.0. R gives positions of extrema. 





R (minima) 




R (maxima) 






3.84 


-0.00002484 


7.02 


-0.00002978 


{a = 0.0 


10.18 


-0.00003270 


13.33 


-0.00003479 


b = 0.1} 


73.1 


-0.00004785 


101.3 


-0.00005033 




98.2 


-0.00005009 


208.1 


-0.00005583 




5.13 


-0.54496245 


8.48 


-0.62000305 


{a = 0.0 


11.74 


-0.66506689 


14.97 


-0.69788474 


b = 0.9} 


18.18 


-0.72371867 


97.3 


-0.94062109 




295.5 


-1.08209885 


198.0 


-1.03113511 
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